Acoustic collective excitations and static dielectric response 
in incommensurate crystals with real order parameter 
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Starting from the basic Landau model for the incommensurate-commensurate materials of the 
class II, we derive the spectrum of collective modes for all (meta)stable states from the corresponding 
phase diagram. It is shown that all incommensurate states posses Goldstone modes with acoustic 
dispersions. The representation in terms of collective modes is also used in the calculation and 
discussion of static dielectric response for systems with the commensurate wave number in the 
center of the Brillouin zone. 
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I. INTRODUCTION 



The uniaxial incommensurate-commensurate (IC) materials of the II class have the wave number of commensurate 
ordering, q c , either at the center (q c = 0) or at the border (q c = ir) of the Brillouin zone. Here the unit length is 
taken equal to the lattice constant. The corresponding Landau free energy functional is then an expansion in terms 
of a real order parameter u(z) 0^], 
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Since the parameter c may acquire negative values, it is necessary to include into the expansion (|l|) the further, second 
derivative term, with a coefficient d presumably positive. 

The mean- field phase diagram follows from the simple variational procedure for the functional It contains 
disordered, commensurate, and various stable and metastable p[ periodic phases. The most important among 

the latter is almost sinusoidal, 
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with weak higher harmonics 

Although represented in terms of the one-component (real) order parameter, the periodic phases like (|2|) are generally 
incommensurate with respect to the underlying lattice. From the other side, it was stated || that incommensurate 
orderings, including those close to the commensurabilities of the II class, should have to be represented by at least 
two quantities, i.e., by the amplitude and the phase of some periodic modulation. This expectation originates from 
the experience with the I class of IC systems, characterized by at least two-dimensional order parameters. More 
specifically, the incommensurate orderings are then most often represented by the modulation of the phase of a 
complex order parameter. Since the model (^) apparently does not comprise a phase variable, it was interpreted 
as a reduction, appropriate only for the description of the commensurate ordering, of some richer physical models 
that include at least two coupled modes, i.e., one-component order parameters. The explicit Landau model with this 
property was formulated by Levanyuk and Sannikov W], and widely explored afterwards fljulpljicjl . The alternative 
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approaches along same lines, attempting primarily to explain the phase diagram for betaine calcium chloride dihydrate 
(BCCD), were also proposed The physical arguments for such general approach to insulators from the II class 

were given by Heine and McConnell |H| . 

Coming back to the model (|l]), it is important to emphasize that, either already bearing the whole physical relevance 
(like in charge density wave materials [fL4|), or being derived from a more complex starting scheme it accounts 
for the cross-over between orderings with real and complex order parameters. In that sense the solutions like (Q) 
are examples of the mean field approximation for the latter. This can be most easily recognized by looking at the 
dispersion of the quadratic part of the expansion (Q) in the reciprocal space, with the biquadratic ("bottle bottom") 
dependence on the wave number (see Eq. ||). This dispersion is an expansion around a pair of symmetry related 
minima, that in addition takes properly into account the symmetry condition on evenness with respect to the center 
(or the boundary) of the Brillouin zone. In this sense the model (Q) can be qualified as the basic one for the II class, 
while the more complex Landau expansions are necessary in, physically possible, but non-generic cases when two 
or more hybridized modes are simultaneously close to instability. 

In order to justify the above statement at the mean field level, one has to prove that incommensurate periodic 
states of the model (|lj) fulfill the crucial general property of incommensurately modulated orderings, namely that for 
each of them there exists a Goldstone mode with the frequency £lo(k) that vanishes at k — 0, and generally has a 
finite slope (phase velocity) dilo(k) / dk = vq in the long wavelength limit k — > 0. Such mode should exist in the whole 
control parameter (e.g. temperature) range of (meta)stability for a given state. It would correspond to standard 
acoustic phason branches for IC orderings in the materials of the I class. The existence of the Goldstone mode in 
incommensurate states like (||) for the materials of the II class is to be contrasted to a widely accepted belief that 
these states do not have a phason mode. 

The requirement tto(k = 0) = guarantees the presence of translational degeneracy of the ordered state. In the 
particular case ([!]) this means that the periodic solutions should be invariant to translations z — > z + zq with arbitrary 
zq. This is obviously fulfilled, since the kernel in the functional ([!]) does not depend explicitly on z. 

It remains to find out whether there exists a Goldstone mode with a nontrivial dispersion [fio(fc) 7^ f° r & 7^ 0], and, 
if so, to determine its properties at the critical lines of the phase diagram for the model ([j]). With this aim, we calculate 
in the present work the spectra of collective modes for all basic orderings that follow from the corresponding Euler- 
Lagrange equation. To this end we generalize the usual eigenvalue problem for collective modes to the systems with 
non-standard free energy densities, in particular to those like (Q) with higher order terms in the gradient expansion. 
Details of this method are presented in Ref. Combining analytic and numerical analysis we show that the 

Goldstone modes for the incommensurate orderings like (^|) exist in the whole range of their stability. Furthermore, 
it appears that collective modes of the model (|l|) have some peculiar, experimentally observable, properties. For 
illustration we discuss the contribution of collective modes to the dielectric response of the II class materials for which 
the order parameter coincides with their electric polarization. 

In Sec. II we continue by considering collective modes for the sinusoidal state (j|) and for the disordered state of the 
free energy (|l|). The expressions for corresponding dielectric susceptibilities are derived in Sec. III. The concluding 
remarks are given in Sec. IV. 



II. COLLECTIVE MODES 



The Euler-Lagrange equation for the model (|l|) reads 

^d 4 u d 2 u 
dz 4 dz 



d— + c—^ + au + bu A = 0. (3) 



Given a solution uq(z) of this equation that also satisfies additional extremalization requirements [ ^6| , the correspond- 
ing eigenvalue problem is defined by 

= dr)n'(z) - ctiq(z) + [a + 3bv%(z)] r)n(z) = n 2 i]n(z). (4) 

The spectrum of collective modes for a given (meta)stable solution uq(z) is defined by those non- negative values of 
il 2 for which the problem (^J) has normalizable solutions t]q(z). For periodic orderings, Uq(z + 2w/Q) — uq(z), we 
adopt the Floquet analysis of this problem |l5| . In particular, it is convenient to use for such orderings the Bloch 
representation T)n(z) — > i] n _k(z) 1 with 

r) n ,k( Z ) = e lkZl $„,k( z ), *n,fe (z + ^] = *n,fe(2). (5) 
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The wave number k is limited to the first Brillouin zone formed by the incommensurate modulation, —Q/2 < k < Q/2, 
and n is the branch index. The formulation (^|) enables the extension of the numerical method, developed for finding 
the solutions uq(z) to the calculation of dispersions f2 n (fc) and Bloch functions ^ n ^k{z). 

While the previous calculations of the spectrum of collective modes for the incommensurate solution u s given by 
Eq. U were approximative |l7| ], the present approach (f|,||) enables exact results, as shown in Fig. 1. The necessary 
condition for the (meta)stability of this and other periodic solutions is c < 0. Then it is convenient to introduce 
the parameter A = ad/c 2 and reduce the model ([j]) to a single parameter problem The lowest among branches 
from Fig. 1, that with the long wavelength dispersion O (fc) = vck, is the Goldstone mode. As is seen from Fig. lc, 
the phase velocity vq tends to zero as A approaches the lower stability edge at \ s = —1.835 (Fig. la), and remains 
finite at the second order phase transition to the disordered state at A^ ~ 0.25 (Fig. lb). In the latter figure we use 
both, Brillouin and extended, schemes for the fc-space, since the former becomes irrelevant for A > Aj(j. The collective 
modes then reduce to the unique dispersion curve for the disordered state, 

il 2 d (k) = a + ck 2 + dk 4 , (6) 

defined by the gradient expansion in Eq. [l]. 

Other collective modes for u s are massive. The dispersion curves for three of them are also shown in Fig. 1. The 
gap of the lowest lying one, f2i(0), tends to zero as A — > Aid- As is seen in Fig. lb, at the very transition to the 
disordered state, A = Ay, this mode has an acoustic dispersion, with the phase velocity equal to vg, i-e., to that of 
the Goldstone mode. 

The spectra of collective modes for other, metastable, periodic solutions of Eq. [| resemble to that of u s . In 
particular, for all of them the phase velocities of Goldstone modes vanish as the parameter A tends towards both 
edges of metastability for a given solution Jlq] , This is to be contrasted to vq for u s , which remains finite at one edge, 
i.e., at A id (Fig. lc). 



III. DIELECTRIC SUSCEPTIBILITY 



The straightforward use of the above Bloch basis (||) leads to the dielectric susceptibility for the incommensurate 
state u s , expressed in terms of branches of collective modes. In the static limit which is under consideration here, it 
is given by 



\9\ 2 



V |^n| 2 
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(7) 



The first term is the response from the Goldstone mode. Here g = jrf^ L (z)dz, and ^^\z) is the coefficient in 
the long wavelength expansion of the corresponding Bloch function (pi), 



* i 1 du s (z) , T (i) , . 
*o,fc(z) =^f7— l^ 1 + k^\z) + ... 
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where J\f is the normalization constant. The expansion of Af in powers of k does not contain the term linear in k 7 i.e., 
in the lowest order approximation it is given by 
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The residual sum in Eq. Q includes the contributions from massive collective modes, i.e., those with finite gaps O n (0). 
Corresponding oscillatory strengths are given by 



n,k=0 



{z)dz\ 



(10) 
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It can be shown |i8| that the coefficient g remains finite as A approaches the stability edge at A s = —1.835, so that 
then cti diverges as Vq . From the other side, this coefficient vanishes together with the function ^i(z) at the transition 
to the disordered state (A = Aj d = 0.25). Since vq remains constant, it becomes clear from the expression (Q) that 
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in this limit the finite contribution to on comes from the residual sum. When calculated directly, by expanding the 
linear response equation in powers of A^ — A, this contribution reads uM 



m~ d 1 



A 



2 (A td - A) 2 

(I-a)(|-a) 



(11) 



The previous calculations HJl!| were limited only to the leading term in this expansion. 

The careful inspection of the residual sum in (f7|) shows that it contributes to the susceptibility ([tl]) only via one 
mode, denoted by in Fig. 1. We note that the oscillatory strength ( [To| ) of the lowest massive mode fii vanishes, so 
that it does not contribute to the susceptibility ([?]). The mode ^2 "survives" the phase transition at Aid, and remains 
the only contribution, presented by the dashed curve in Fig. lb. The dielectric susceptibility in the disordered state 
at A > \id is given by 

f^(0) a 

Thus, at the phase transition from the incommensurate to the disordered state the dielectric susceptibility varies 
continuously, but with a finite jump in at A = Aid. 

To summarize the above discussion of the dielectric responses on and ay, we show in Fig. 2 their dependence on 
the parameter a (i.e. on temperature). Other parameters are fixed. Fig. 2 also includes the susceptibility for the 
commensurate solution u c = ±^/—a/b, which is stable in the ranges A < — | for c < 0, and a < for c > 0. It 
is separated from the disordered state Ud by the second order transition at the line a — 0, c > 0, and from the 
incommensurate state u s by the first order transition, defined by the line c < 0, A —1.177. 

The static susceptibility of the commensurate state u c is given by 

111 

(13) 



a + 3bu 2 c Q2( ) 2a ' 

where 

fl 2 c (k) = dk 4 + ck 2 -2a (14) 
is the dispersion for the corresponding unique collective mode that follows from Eq. H after inserting uq(z) — u c . We 



see from Eqs. (12) and ( |13| ) that the static susceptibility shows a standard type of divergence at the line of the second 
order transition from the disordered to the commensurate state, a — 0, c > 0. From the other side, a c has a finite 
value at the stability edge (c < 0, A c = — |) for u c in the regime of coexistence with the incommensurate state u s . 
Thus, we propose an asymmetric behavior of the susceptibility as one passes through the hysteresis range (in, e.g., 
temperature), bounded by the values A c and X s in the parameter A. Namely, as shown in Fig. 2, by cooling through 
the incommensurate state one ends with the divergence of on before passing into the commensurate state at A s , while 
by heating through the commensurate state one passes into the incommensurate state at A c with the jump in the 
susceptibility equal to [a c — ai]\ = \ c . 

Due to the existence of other metastable periodic solutions [Q in the above coexistence range of the parameter A, 
the temperature variation of the susceptibility may be even more complicated than that schematically presented in 
Fig. 2. The more detailed analysis ]l5[ | shows that, provided that for a given material and in particular circumstances 
some of these states are stabilized, the static susceptibility will diverge at both edges of their stability ranges, again, 
like in the case of the state u s , due to the vanishing of phase velocity for the corresponding Goldstone modes. Since 
these stability ranges are rather narrow in comparison with A c — A s ~ 1.71, the stabilization of highly nonsinusoidal 
periodic states fl would be manifested by dramatic variation of susceptibility in relatively short temperature intervals. 



IV. CONCLUSION 

In the above analysis it was clearly shown that the model (Q), although formulated in terms of a single and real 
order parameter, provides the existence of the acoustic Goldstone mode for the incommensurate ordering. For small 
wave numbers (k <C Q) this mode simply generates compressions and dilatations of the sinusoidal modulation u s (z), 
i.e., its physical content is same as that of standard phason mode for incommensurate orderings in the II class of 
IC systems. Amplitude fluctuations of this modulation are generated through higher, massive modes from Fig. 1. 
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Still, the lowest among these massive modes have some peculiar properties. For one of them the gap tends to zero by 
approaching the transition from, e.g., the sinusoidal state to the disordered state. However this mode is not optically 
active. The whole optical activity at this transition comes from the another massive mode in order, as shown in Fig. 1. 

In contrast to standard approaches |2(|, the dielectric functions for all (meta)stable states are here strictly repre- 
sented in terms of collective modes. We are thus able to distinguish modes with finite contributions in the optical 
response from those with no dipolar polarization. Note that the present analysis of dielectric response is limited to 
systems with uniform commensurate orderings (those with q c — 0). The extension to systems with dimcrizations, as 
well as the equivalent treatment of collective modes and dielectric response of extended Landau models for IC systems 
of class I pH , are under current investigations. 

The above analysis shows that incommensurate states that follow already from the simplest basic version of the 
model (|l|) for the IC systems of class II, have the same physical properties as those of class I, represented in terms 
of the complex order parameter. This model is still insufficient for the explanation of rich phase diagrams of some 
well-known representatives of class II, like thiourea and BCCD. In this respect the question which remains is, how 
to make an appropriate extension without invoking an additional order parameter, in a way analogous to the recent 
proposal for the class I Q , that would as well stabilize other commensurate states with wave numbers close to q c = 
or q c = 7r / 'a. 
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FIG. 1. 
lines), A = 
(Fig. lb). 



The dispersion curves for the lowest collective modes of the incommensurate state (g), for A = \ s — — 1.8 (full 
— 0.5 (dashed lines) (Fig. la), and for A = A^ = 0.25 in the Brillouin (full lines) and extended (dashed line) schemes 
The dependence of the phase velocity of the Goldstone mode on the parameter A is shown in Fig. lc. 




G 



1.5 



1.0 



C5 



0.5 



0.0 
-2.0 



-1.0 



0.0 



1.0 



FIG. 2. Dielectric susceptibilities for the incommensurate state (full line), commensurate state (dashed line), and disordered 
state (dotted line) as functions of the parameter A. The critical values of A shown in the figure, are introduced in the text. 
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